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Human-wildlife conflict is a leading cause of adult mortality for large carnivores worldwide. Train collision 
is the primary cause of mortality for threatened grizzly bears (Ursus arctos) in Banff National Park. We 
investigated the use of stable isotope analysis as a tool for identifying bears that use the railway in Banff. 
Rail-associated bears had higher (5^^N and values than bears sampled away from the rail, but similar 
(5^^C values. Because elevated (5^^N values are indicative of higher animal protein consumption, rail- 
associated bears likely preyed on ungulates that foraged along the rail or scavenged on train-killed animals. 
The higher (5^^S values in bear hair could have resulted from bears consuming sulfur pellets spilled on the 
rail or through the uptake of sulfur in the plants bears or animals consumed. Similar S^^C values suggest 
that the two types of bears had generally similar plant-based diets. Results from this study suggest that 
stable isotopes analysis could be used as a non-invasive, affordable, and efficient technique to identify and 
monitor bears that forage on the railway in Banff and potentially other transportation corridors worldwide. 
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1. Introduction 

Large-bodied carnivores typically range over large areas [1], exhibit slow population growth rates 
[2], and occur at low population densities [3]. All of these factors increase their susceptibility 
to extirpation, particularly in areas with increasing anthropogenic impacts [4]. Consequently, 
human-carnivore conflict is often the leading cause of adult mortality for large carnivores, 
even in protected areas [5]. Preventing or mitigating human-wildlife conflict is dependent on 
a mechanistic and behavior-based understanding of the human-wildlife interface [6]. In partic- 
ular, it is important to identify the individuals within a population that exhibit behavior(s) that 
increase their risk of human-induced mortality and target the associated individuals and specific 
behavior(s) that engender conflict with the appropriate monitoring and management practices. 
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Figure 1 . The transportation corridor in Banff National Park, Alberta, Canada. 

Traditional methods for monitoring such animals through time generally rely on capture and 
telemetry, but these techniques can be invasive, costly, and time-consuming [7]. Meanwhile, 
the rate of human-wildlife conflict is increasing worldwide [8], creating an imperative to iden- 
tify novel methods that are non-invasive, inexpensive, and efficient at identifying conflict-prone 
individuals. 

Human-carnivore conflict occurs frequently along the transportation corridor in Banff National 
Park (hereafter, Banff) [9]. Carnivores are killed in Banff by automobiles on the Trans-Canada 
Highway (TCH) [10,11] and by trains on the Canadian Pacific Railway (CPR; Figure 1) [12]. 
Grizzly bears (Ursus arctos) have attracted particular concern because the local population repro- 
duces slowly and has one of the lowest densities in North America [13]. High rates of mortality 
caused by collisions with automobiles on the TCH have been reduced by the installation of exclu- 
sionary fences [10] and wildlife crossing structures [9]. Unfortunately, the rate of collisions with 
trains on the CPR is now the leading cause of grizzly bear mortality in Banff [14]. No one knows 
what change in bear or train behavior has caused this increase, but identifying conflict-prone 
individuals could make it possible to mitigate this source of mortality. 

The most prevalent hypotheses for the increase in train-caused mortality of grizzly bears in Banff 
National Park apply to transportation systems everywhere. Many species of animals are attracted 
to roads for the foraging and travel opportunities they provide, while lacking the evolutionary 
heritage to identify vehicles as mortality threats. The high rate of attraction to roads across 
taxonomic groups contributes to the role of mortality as potentially being the most detrimental 
impact of all road effects [15-17]. Although there has been much less research on the effects of 
railways on wild animal population, it is clear that grizzly bears in Banff use the rail both for 
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travel and food [12,18]. Grizzly bears that use the railway in Banff forage on plants and a variety 
of human-derived foods including train-killed animals and deposits of agricultural products being 
carried by the rail (hereafter, grains) [12,18,19], including wheat (Triticum aestivum), barley 
(Hordeum vulgare), canola (Brassica spp.), lentils (Lens culinaris), peas (Pisum sativum), and 
com (Zea mays) [12,20]. In addition to grains, researchers and park staff frequently observe sulfur 
pellets along the tracks and occasionally in bear scats (J. Hopkins, personal observation). Grains 
and sulfur pellets are deposited on the rail during train derailments, leaked from grain hopper cars, 
or blown off cars after being spilled on their exteriors during loading [19-22]. Not all grizzlies 
in the Mountain Parks appear to forage on the rail (Hal Morrison, personal communication), but 
those that do are presumed to be more susceptible to train-caused mortality [21,22]. 

Grizzly bears that are at risk of being killed by trains on the CPR in Banff epitomize the kind 
of human-wildlife conflict for which it would be helpful to have a non-invasive, inexpensive, 
and efficient method for identifying conflict-prone individuals. The goal of this study was to 
determine if stable isotope ratios could be used for such a purpose. Stable isotopes ratios have 
been used in past studies to identify individuals that forage on human-derived foods [23-25]. For 
instance, black bears (Ursus americanus) that foraged on meat-rich, human foods [25] or animals 
[26] exhibited higher nitrogen isotope ratios (^^N/^'^N, expressed as 5^^N values) than bears that 
primarily foraged on plants. Similarly, individuals that foraged on C4 plants such as corn from 
agricultural fields or corn-derived products, which include processed foods or corn-fed livestock 
[27,28], had higher carbon isotope ratios (^^C/^^C, expressed as 8^^C values) than individuals 
that forage on C3 plants exclusively [23,25]. 

We hypothesized that the isotopic composition of hair for bears killed on or captured near the rail 
(hereafter, rail-associated bears) in Banff would reflect the anthropogenic foods they likely forage 
on. Our hypothesis predicted: (1) higher 5^^N values for rail- associated bears due to preying on 
ungulates that foraged along the rail or scavenging on rail-killed animals; (2) higher ^^'^S values 
for rail-associated bears due to ingesting sulfur pellets directly (e.g. when lapping up grains) 
or indirectly (e.g. via plants that grew along the rail in sulfur-enriched soil or via animals that 
consumed those plants or pellets); and (3) different 5^^C values for rail-associated bears due to 
foraging on different plants (e.g. growing along the ballast) or train-spilled agricultural products 
(e.g. com elevating 8^^C in their tissues). We tested these predictions by comparing the isotope 
values in hair of rail- associated bears to those of bears sampled away from the rail throughout 
the park. 



2. Materials and methods 
2.1. Study area 

Banff National Park (51.2°N, 115.5°W; 6,641 km^) is located approximately 150km west of 
Calgary, Alberta (Figure 1) in the Canadian Rocky Mountains. Banff receives more than three 
million visitors each year. The park contains rugged mountains with elevations ranging from 
1000 to 3500 m with a climate characterized by long, cold winters, and short, cool summers [29]. 
Montane and subalpine ecoregions in Banff are dominated by lodgepole pine (Pinus contorta), 
Douglas-fir (Pseudotsuga menziesii), white spruce (Picea glauca), Englemann spruce (Picea 
englemannii), trembling aspen (Populus tremuloides), and natural grasslands. Primary grizzly 
bear foods include sweetvetch roots (Hedysarum spp.), bearberries (Arctostaphylos uva-ursi), 
graminoids (Family), horsetails (Equisetum spp.), buffalo berries (Shepherdia Canadensis), and 
other fruits (Vaccinium spp.) [30]. 

In addition to grizzly bears, Banff hosts a diversity of medium- to large-bodied mammals 
including black bears, wolves (Canis lupus), coyotes (Canis latrans), foxes (Vulpes vulpes). 
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cougars {Puma concolor), lynx (Lynx canadensis), wolverines (Gulo gulo), moose (Alces alces), 
elk (Cervus canadensis), mule deer (Odocoileus hemionus), white-tailed deer (O. virginianus), 
bighorn sheep (Ovis canadensis), and mountain goats (Oreamnos americanus). Most of these 
species have been killed on the TCH and the CPR (~82 km) over the past 20 years in Banff [12]. 

2.2. Sampling 

We used isotope values from three different samples of bear hair collected in Banff in 2000-2012. 
We collected hair samples from independent grizzly bears killed on the tracks in 2000-2011 
(n = 9) and those captured adjacent to the railway in 2012 (/i = 9) (Table 1). We also used hair 
samples that were collected away from the rail at wildlife crossing structures, hair traps, and bear 
rubs in 2006-2008 (see Sawaya et al. [31] for detailed descriptions of sampling methods). 

Parks Canada personnel collected hair from grizzly bears killed on and captured near the rail. 
Park personnel visited each mortality site and collected an envelope of hair from each dead bear. 
Some personnel noted that the stomachs of train-killed bears contained grain (Parks Canada, 



Table 1. Isotope values {S^ C, (5^N, (5^^S) for bears killed and captured on the railway in Banff 
National Park, 2000-2012. 



Park ID 


Sex 


Date killed 


Hair-year 




(^i^N 


8^^S 


Killed bears {n = 9) 














171 


F 


8/17/00 


1999 


-23.7 


4.5 


10.9 


46 


F 


9/13/04 


2004 


-22.6 


3.8 


6.9 


4202^ 


UNK 


5/31/05 


2004 


-23.5 


4.6 


9.8 


66^ 


F 


8/19/05 


2003 


-23.2 


6.3 


6.0 


4346 


F 


6/18/07 


2006 


-23.2 


4.3 


5.6 


109 


F 


5/27/10 


2009 


-22.7 


4.1 


6.8 


4586 


M 


6/25/10 


2009 


-22.5 


5.1 


7.2 


176 


F 


9/30/11 


2001 


-23.6 


4.7 


10.8 


UNK 


M 


11/24/11 


2011 


-22.6 


5.7 


8.6 


Mean 








-23.1 


4.8 


8.1 


SD 








0.5 


0.8 


2.0 


Captured bears (n = 9) 














122 


M 


4/18/12 


2011 


-23.0 


7.0 


9.0 


125 


M 


5/23/12 


2011 


-23.1 


3.3 


7.5 


126 


M 


5/26/12 


2011 


-23.7 


5.4 


9.2 


72 


F 


6/1/12 


2011 


-23.7 


5.4 


10.2 


130 


F 


6/7/12 


2011 


-23.6 


4.7 


3.2 


132 


M 


6/14/12 


2011 


-23.4 


5.0 


4.6 


134 


M 


6/23/12 


2011 


-23.1 


6.0 


6.8 


135 


F 


6/25/12 


2011 


-23.5 


4.7 


2.4 


136 


M 


6/27/12 


2011 


-23.3 


6.7 


5.7 


Mean 








-23.4 


5.4 


6.5 


SD 








0.3 


1.1 


2.7 


Males (n = 8) 














Mean 




-23.1 


5.5 


7.3 






SD 




0.4 


1.2 


1.6 






Females(w = 9) 














Mean 




-23.3 


4.7 


7.0 






SD 




0.4 


0.7 


3.1 






Rail-associated bears (n = 18) 














Mean 








-23.2 


5.1 


7.3 


SD 








0.4 


1.0 


2.5 



Notes: All isotope values are expressed in delta (8) notation as parts per thousand (%o). Hair-year denotes the year 
isotopes were assimilated into the hair. Suess corrections account for the global decrease of ^^C concentration in 
Earth surface carbon reservoirs, primarily due to fossil fuel burning over the past 150 years [42]. 
^Hair collected a year prior to being killed. 
Sub-adult bear; age and sex unknown. 



326 



J.B. Hopkins 111 et al. 



unpublished data). Parks Canada personnel also collected hair from immobilized bears captured 
near the railway. Park personnel captured grizzly bears using a combination of free-range darting 
and culvert traps using protocols approved by Parks Canada Animal Care Committee and in 
accordance with guidelines approved by the American Society of Mammalogists [32]. The capture 
team immobilized bears using low impact, slow-inject darts containing 3 mg/kg of tiletamine- 
zolazepam and 2 mg/kg of xylazine. Once immobilized, the capture team recorded standard 
morphological measurements for each bear, fit it with a GPS radio-collar, injected 0.2 mg/kg 
of atipamizol to reverse the effects of xylazine, and then released the bear. We sub-sampled >2 
guard hairs of equal length from each of these samples. We included only those bears in our 
samples that had been independent of their mothers for more than one year as nursing may enrich 
^^N and deplete ^^C in the tissues of young bears [33] (Table 1). Because grizzly bears molt 
during spring and summer months, we assumed full-length guard hairs collected from captured 
bears in spring 2012 were grown and contained isotopes derived from bear foods consumed 
in 2011 [34,35]. 

Recent studies collected hair samples in 2006-2008 to estimate the demographic rates of 
grizzly bears in Banff [31] and to evaluate population connectivity using genetic data [36,37]. We 
identified bears for isotope analysis = 1 15) by querying unpublished data from Sawaya et al. 
[31]. First, we sorted individuals by year and by hair availability. Next, we queried samples by 
season. We used only guard hairs that were of full and equal length and included only those samples 
collected in May- August because samples collected in September onward were likely grown that 
year. We considered individuals that were recaptured in successive years to be independent samples 
to acknowledge potential differences in diet among years. 



2.3. Sample preparation and stable isotope analysis 

We rinsed hairs with 2:1 chloroform-methanol solution to remove surface oils and then air-dried 
the samples for stable isotope analysis (SI A). We conducted SI A at University of Alberta (hereafter, 
U Alberta Lab) and University of Ottawa (hereafter, U Ottawa Lab). The U Alberta Lab weighed 
single hairs into tin capsules (4 mmx 6 mm; #D1006 from Elemental Microanalysis, Okehampton, 
UK) using a Starorius CP2P micro-analytical balance. The U Alberta Lab analyzed samples for 
their carbon (^^C/^^C) and nitrogen (^^N/^^N) isotope compositions by continuous flow methods 
using a Costech 4010 HAS Elemental Analyzer (Costech Analytical Technologies Inc., Valencia, 
USA) coupled to a ThermoFinnigan Conflo III and ThermoFinnigan Delta+Advantage stable 
isotope mass spectrometer (Thermo-Scientific, Waltham, USA). The U Ottawa Lab analyzed 
samples for their sulfur (^^S/^^S) isotope compositions by loading the encapsulated hairs into a 
Vario Micro Cube (Elementar Americas, Mt. Laurel, USA), coupled to a DeltaXP IRMS (Thermo- 
Scientific). 

Stable isotope ratios are expressed in delta (5) notation as parts per thousand or per mil (%o). 
SIA results are reported and defined as 

_ ('^^/^^) sample ^ 
(^X/^X) standard 

where is the heavier isotope (^^C, ^^N, ^"^S), and the lighter isotope (^^C, ^"^N, ^^S) in 
the sample (numerator) and international measurement standard (denominator) [38]. Reference 
standards are V-PDB for carbon (ib 0.2 %o), atmospheric N2 for nitrogen (±0.2 %o), and Vienna 
Canon Diablo Troilite for ^^^S (ib 0.2 %o) [39-41]. 

We corrected 5^^C values for the Suess effect, the global decrease in ^^C in Earth's atmospheric 
CO2 due to fossil fuel burning over the past 150 years [42]. Following Chamberlain et al. [43], 
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we applied a time-dependent correction of -0.022 %o per year through 2011 for all hair samples 
except those collected in 2012. 

2.4. Statistical analyzes 

We calculated intra-sample error (x, ISD) by conducting isotope analysis on multiple hairs of 
equal length within samples. This error (^^^N: x = 0.3 ib 0.2%o; ^^^S; x = 1.0 =b 1.0%o; ^^^C: 
X = 0.2 lb 0.2%o; Table SI) accounts for the inherent isotopic variation among hairs within a 
sample and acknowledges the degree of precision afforded by conducting SIA (Table SI). We 
conducted our analyzes described below using mean isotope values for a minimum of two hairs 
from each killed and captured bear and one hair from each hair-snared bear. 

We tested for differences in isotope values between bears killed on and captured near the rail, 
between male and female bears, and between rail-associated bears and those sampled away from 
the rail throughout the park. We used r-tests to assess for differences in ^^'^S because these isotope 
values follow a normal distribution. We also used Mann-Whitney, Wilcoxon tests to test for 
differences in both 3^^C and 5^^N because these isotope values follow non-normal distributions. 
We used R (2.15.1) to conduct all significance tests and set a = 0.05. 



3. Results 

We conducted univariate analyses to identify differences in isotope values for the two samples of 
rail-associated bears and for the two sexes in both types of samples. As expected, we found that the 
18 grizzly bears killed on or captured near the railway had similar isotope values ((5^^N: W = 24, 
P = 0.16; (5^4s: t = 1.37, P = 0.19; ^^^C: W = 55,P = 0.21; Table 1), suggesting they had sim- 
ilar diets. Likewise, isotope values were similar for male and female bears for rail-associated 
samples (^^^N: W = 55.5, P = 0.07; ^^^S; t = 0.29, P = 0.78; ^^^C: W = 49.5, P = 0.21; 
Table 1) and hair-snared samples (^^^N^aie = 4.1, ^^^Nfemaie = 3.7, W = 1423.5, P = 0.256; 
8^^Sm.ic = 5.5, 8^%,m,ic = 5.2, t = -0.63, P = 0.53; 8'^C^,i, = -23.2, ^^^Cfemaie = -23.2, 
W = 1704.5, P = 0.65; Table S2). Accordingly, we aggregated rail-associated bears into a single 
group. 

When we compared the isotope values for rail-associated bears and those sampled away from 
the rail, we found significant differences for both 8^^N and 8^^S values (Figure 2(a) and 2(b)) but no 
difference in 5 ^^C values (Figure 2(c)). Mean values of 5 ^^N were 1 .2 %o greater for rail-associated 
bears relative to hair-snared bears (Tables 1 and S2; W = 1700.0, P =< 0.001). Likewise, the 
values of ^^'^S were 1.9 %o greater for rail-associated bears (Tables 1 and S2; W = 1490.0, P = 
0.003). By contrast, mean values of 8^^C were effectively equal (W = 992, P = 0.78) for the two 
groups of bears (Tables 1 and S2). 



4. Discussion 

Rail-associated mortality of grizzly bears has increased dramatically in Banff National Park. 
Evidence from this study suggests that it is plausible that SIA of hair can be used as a reliable 
tool to reveal the rail-based foraging behavior of grizzly bears in Banff. Our results demonstrated 
that rail-associated bears had different isotope values than bears sampled away from the rail, 
supporting the hypothesis that hair samples could be used to identify individuals that forage on 
the rail and are thereby at risk of being killed by trains. Rail-associated bears in our study had higher 
8^^N and ^^'^S values than hair-snared bears, but similar 8^^C values. We found no differences 
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Figure 2. Boxplots depicting nitrogen ((5^^N; Panel a), sulfur (^^'^S; Panel b), and carbon ((5^^C; Panel c) isotope 
values derived from the hair of grizzly bears killed (black squares) or captured (black circles) on the railway (defined as 
rail-associated bears; Panel d) and those sampled via hair-snare (grey diamonds; Panel d) in Banff National Park, Alberta, 
Canada, 2000-2012. 



in isotope values between bears that were killed on or captured near the rail or between the two 
sexes, suggesting that the nitrogen and sulfur isotopes contained in hair are robust measures for 
detecting independent grizzly bears that forage on the rail in Banff. 

Higher (5 ^^N values measured in the hair of rail- associated bears are indicative of diets that were 
higher in animal protein [26,44,45]. Grizzly bears in Banff are considered to be protein-limited 
[13]; although they periodically feed on ungulates, ants are a more consistent source of animal 
protein [30]. It is plausible that rail- associated bears had access to more ants along the rail in the 
Bow Valley, but it is more likely that higher animal protein in the diets of rail-associated bears was 
due to preying on ungulates or scavenging on rail-killed animals. Dozens of animals are killed on 
the rail each year in Banff (Parks Canada, unpublished data), and many large carcasses are found 
and removed by park staff. Transplanted carcasses are deposited in nearby restricted areas where 
they are available to carnivores, and many bears that use the rail are also known to feed at these 
sites (D. Garrow, Parks Canada, personal communication). Whether feeding on carrion on the 
railway or at carcass dumps, bears that used the transportation corridor in Banff likely obtained 
substantial protein subsidies by feeding on animals killed by trains, thus elevating 5^^N values 
in their hair [26]. Even though (5^^N values for rail- associated bears are significantly higher than 
hair-snared bears, the actual difference between groups might be greater if we included dependent 
bears that nursed during the previous year in our rail- associated sample or excluded these bears 
from our sample of hair-snared bears [33]. 
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Rail-associated bears also had higher ^^'^S values than hair-snared bears, indicating that some 
bears may have fed directly on sulfur pellets or foraged on plants or animals that assimilated ^"^S 
into their tissues. Sulfur isotopes may prove to be a particularly informative predictor of railway 
use by bears and other animals in the area because ^"^S occurs in relatively low abundance in 
bear foods; one exception occurs in whitebark pine seeds [35], but these are not common in the 
mountain parks (J. Whittington, personal observation). 

Although it appears that ^^'^S is a clear indicator of rail- associated foraging for grizzly bears 
in Banff (Figure 2(b)), it is important to determine the source(s) of ^"^S in Banff and the foraging 
behavior(s) bears use to acquire it. As reflected by the historic attraction of railway passengers to 
Banff afforded by the Sulfur Mountain hot springs, ^"^S may occur at elevated levels naturally in and 
around the rail corridor. In addition, soil adjacent to the rail was augmented with ^"^S for decades 
when railway firemen shoveled the spent coal from steam engines onto the ground adjacent to 
the rail (K. Roberge, personal communication). Modem combustion engines fueled with diesel or 
gasoline emit sulfur dioxide, which could also result in heightened 5^^S values in the vegetation 
[46] parallel to the TCH and CPR in Banff. Whether from geological or anthropogenic origins of 
historic or recent eras, plants growing in those areas may have elevated 5^^S values, which might 
then be consumed by bears directly or indirectly through the animals they feed on. Finally, some 
observers in Banff have detected sulfur pellets in the scats of bears, which may be consumed 
intentionally or incidentally when bears target other products on the rail. Much further work will 
be needed to estimate the proportion of sulfur isotopes emanating from these potential sources. 

The lack of difference in the carbon isotopes of rail-associated and hair-snared bears 
(Figure 2(c)) was contrary to our prediction and surprising because corn, a C4 plant that is rela- 
tively high in ^^C, is transported on the railway. Moreover, some rail-associated bears in Banff 
have been known to consume other sources of anthropogenic foods, which have resulted in higher 
values of hair for black bears elsewhere [23,47]. We may have failed to detect a difference 
because none exists, which would be expected if com were rarely deposited on the tracks, or if 
a tendency to consume other forms of anthropogenic food is uncommon in the rail-associated 
bears. Without these influences, similar 5^^C values would be expected as grizzly bears in Banff 
have homogeneous plant-based diets [30]. Although it was beyond the scope of this study, inte- 
grating carbon isotopes with other isotopes could refine the differences in plant-based diets in 
rail-associated bears ([26,45,47]; Figure 2(d)). 

Further work that integrates all pair- wise combinations of the three isotope measures and behav- 
ioral data may reveal additional information about the diets of individuals that associate with the 
rail in Banff and other transportation systems. In other contexts, the combined signals of multiple 
isotope values have revealed intra-population niche variation [48-50] and dietary specialization 
of anthropogenic foods [23-25]. Evaluating these differences more comprehensively may make 
it possible to distinguish between rail- associated bears that forage heavily on grains and those 
that forage on other foods while traveling along the rail. In subsequent analyses, greater differ- 
entiation of rail-associated bears that were known to specialize on grain may reveal 5^^C to be 
an informative predictor of this behavior (as suggested in Figure 2(d)), in turn supporting more 
specific management intervention. 



5. Conclusion 

Results from this study demonstrate that stable isotopes in the hair of grizzly bears can distin- 
guish between individual bears that were known to be associated with the railway in Banff and 
those without a known association. Because stable isotopes are ubiquitous in the environment 
[51] and their analyses are standardized [52], this result has important implications for wildlife 
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conservation in other jurisdictions around the world where train colHsions and other transporta- 
tion effects threaten populations. The non-invasive and inexpensive attributes of this method 
make it particularly well suited for the developing nations where rail-based transportation is 
prevalent, but lack of funding precludes more expensive methods for identifying at-risk behaviors, 
populations, and individuals. Once vulnerable individuals are identified, managers could analyze 
their hair repeatedly through time to monitor their behavior and to determine the success of 
mitigation efforts, adjusting their interventions subsequently within a framework of adaptive 
management. Such a behavior-based approach to management could increase the proportions 
over time of both behaviors and individuals that coexist most sustainably with the burgeoning 
human population. 
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